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Ordered phases in Landau paradigm can be diagnosed by a local order parameter, whereas topo¬ 
logically ordered phases cannot be detected in such a way. In this paper, we propose long-range 
mutual information(LRMI) as a unified diagnostic for both conventional long-range order and topo¬ 
logical order. Using the LRMI, we characterize orders in n -|- ID gapped systems as m-membrane 
condensates with 0 < m < n — 1. The familiar conventional order and 2-l-lD topological orders 
are respectively identified as 0-membrane and 1-membrane condensates. We propose and study 
the topological uncertainty principle, which describes the non-commuting nature of non-local order 
parameters in topological orders. 


Introduction - In the Landau paradigm, conventional 
orders (CO) are characterized by the phenomenon of 
spontaneous symmetry breaking. In the thermodynamic 
limit, some symmetry-breaking local operator-the order 
parameter- attains an expectation value. Alternatively, 
one can also consider the symmetry-preserving ground 
state and describe CO by a nonvanishing long-range cor¬ 
relation between order parameters, as was proposed in 
the case of off-diagonal long-range order [1]. Another kind 
of order that is beyond the Landau paradigm, the topo¬ 
logical order (TO), has been studied since the discovery 
of fractional quantum Hall effect [5]. TO’s are usually 
characterized by different theoretical frameworks such as 
topological ground state degeneracy or braiding statis¬ 
tics of the topological quasi-particles. While non-local 
order parameters have been proposed in some particular 
TO’s such as Laughlin state[3] and Z 2 spin liquid [HIS], 
much less is known about such order parameters in gen¬ 
eral TO’s. Since systems with CO and TO both have 
robust ground state degeneracy in the thermodynamic 
limit, it is natural to ask whether they can be described 
in a unified framework. Such a framework is likely to pro¬ 
vide new insight on the definition and characterization of 
TO in higher dimensions, which is much less understood 
than that in two spatial dimensions. 

In this letter, we propose a new measure of long- 
range entanglement, the long-range mutual information 
(LRMI), as a unified language to describe both CO and 
TO. Mutual information between two regions of a many- 
body system is known to upper bound all connected cor¬ 
relation functions them[B]. Therefore, the LRMI between 
far away regions with different topology should be a good 
indicator of both CO and TO. In the following, we will 
provide the definition of the LRMI. Then we shall study 
the LRMI diagnostic of the CO and different TOs in 
2-l-lD, as well as in higher dimensions. As it shall be¬ 
come clear, we view the CO as a condensate of particles, 
2-l-lD TO as a condensate of strings, and higher dimen¬ 
sional TO as that of extended objects of various dimen¬ 
sions. Therefore, the topology of the chosen regions plays 


an important role in understanding the nature of the un¬ 
derlying order. The LRMI description also explains the 
key difference between TO and CO. The non-local or¬ 
der parameters in TO can be defined on intersecting re¬ 
gions, which generically do not commute with each other. 
The non-commuting nature is characterized by a “topo¬ 
logical uncertainty principle”, which states that a pair 
of LRMI on these intersecting regions cannot vanish si¬ 
multaneously. The topological uncertainty principle is 
a smoking-gun indication of long-range quantum entan¬ 
glement in TO. While we will mainly focus on gapped 
systems with a finite correlation length ^ in this letter, 
we shall also provide some comments on gapless systems 
as well. 

Definition of the LRMI - For a given quantum 
many-body state, we denote pA as the reduced density 
matrix of a region A. The entanglement entropy of A is 
given by Sa = S{pa) = —TrpyilogpA- For two disjoint 
regions A and B, their mutual information is defined as 

Ia,b = Sa + Sb — Saub- (1) 

For our purpose, we always take regions A and B to have 
the same topology (denoted as X) and take the limit 
in which the distance between A and B{dist{A, B)) is 
infinite. In this limit, we denote the asymptotic value of 
Ia,b as I{X), and refer to it as the LRMI [5^. There are 
different choices of X that characterize different types of 
orders. For example, on a two-dimensional torus, X can 
be a disk or a cylinder. As we shall see later, these choices 
work as a probe to detect CO and TO respectively. 

LRMI in conventional order - To determine the 
CO in a n -I- ID system, we only need to consider regions 
A and B that are disks I?" with size r ^ The LRMI 
I{D^) is defined in the limit dist(A, R) —>■ 00 with r 
fixed. Given the ground state Hilbert space TLg of a 
gapped system, we propose: 

A gapped system possesses CO when the LRMI 

1(1?") does not vanish for generic states in Tia- 
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In disordered phases, all correlation functions factorize 
at large distance, and so does the reduced density ma¬ 
trix paub- As such, LRMI should vanish. Our justifica¬ 
tion for this proposal is based on the following argument. 
Consider a basis given by the classical states 
with broken symmetries and no connected long-range 
correlation. For example, in a simple Ising ferromagnet, 
the ground state Hilbert space is two-dimensional, and 
the two basis states have all spins up or all spins down. 
A generic ground state is a superposition {-ip) = ipi\i). 
Since the classical states are macroscopically distinguish¬ 
able, the reduced density matrix of \tp) on any small re¬ 
gion X compared to the system size is block diagonal: 

/ Ibil^pi.x \ 

px = I lb 2 pP 2 ,x I Here pi^x is the reduced 


distinguished by non-local processes, e.g., by winding a 
quasiparticle around the torus. Therefore, it is natural 
to expect nontrivial LRMI between non-contractible re¬ 
gions, which characterizes string-like order parameters 
wrapping around the non-contractible loops. In the fol¬ 
lowing, we focus on the ground state Hilbert space T-La 
on the torus of size x Ly. The topology of the 
regions A and B are taken to be x where is an 
winding around a non-contractible prime 1 -cycle 7 on 
the X^ (Here a cycle 7 is defined to be prime when there 
exists no other cycle such that a multiple g (q > 1 ) of it 
is topologically equivalent to 7 ). The linear size r of the 
part is fixed and satisfies ^ <C r <C The LMRI 
I{D^ xTp^) is defined in the limit where dist(H, H) —)■ 00 
with d\si{A, B)/ fixed. We propose 


density matrix of the classical state |i) on the region X. 
Taking X to be the regions H, B and AUH in the LRMI, 
we obtain X(X") = log |^,|2, 

where Ii{D^) is the LRMI in the classical state |i). Fur¬ 
thermore, Xi(X)") = 0 since all the connected correlation 
functions of the classical state |i) vanish at large dis¬ 
tances. Therefore, we obtain 


X(X”) = -^|V>.plog|V-.P, (2) 

i 


which is generically non-vanishing in agreement with our 
proposal. Moreover, it is worth mentioning that the mu¬ 
tual information Ia.b serves as an upper bound for all 
connected correlation functions [S]: 


\{OaOb)^ 

^^""- 2110 ^ 1121105112 ’ 


(3) 


for any operators Oa and Ob supported on regions A and 
B respectively. Here, the subscript “c” stands for con¬ 
nected correlation function and || • || for operator 1 -norm. 
Hence, any CO that is signified by the non-vanishing 
long-range connected correlation function of the order 
parameters has a finite value of X(X)"). 

It is worth noting that X(X>"’) in fact takes the form 
of the Shannon entropy with the probability distribution 
la particular, it can achieve a minimal 
value of 0 but only on classical states. This property can 
be in turn viewed as the definition of the classical states 
in our LRMI approach, which will later connect to the 
discussion of the TO. For systems with spontaneous bro¬ 
ken continuous symmetries, dim('HG) is infinite. While 
the non-vanishing LRMI for generic states in Hq is still 
guaranteed by Eq. its specific form is left for future 
work. 

LRMI in 2+lD topological order - In a general 
2 -I- ID system with TO, all the connected correlation 
functions of local operators vanish at large distance, and 
so does the LRMI between disks(X(X)^)). However, TO 
on a torus has degenerate ground states, which can be 


A gapped 2+lD system possesses TO (but not CO) 
when the LRMI x T^) is finite for generic 

states in Be with a given non-contractible prime 
1-cycle^ andT{D^) vanishes throughout Bg ■ 

We first consider the thin-torus with LyjL^ 1 limit, 
in which our proposal can be succinctly explained. The 
TO on a thin torus can be viewed as the CO of a quasi- 
(1-|-1)D system. The most well-known example is the 
Laughlin fractional quantum Hall state, for which the 
different topologically degenerate ground states on the 
X^ cross over to different charge density wave patterns 
on the thin torus [71-fT^. Another example is given by the 
2-l-lD toric code model which, in the thin-cylinder limit, 
reduces to a spin chain with Z 2 x Z 2 global symmetry. 
As a quasi 1 -l-lD system, this spin chain exhibits CO 
with the order parameter given by Wilson loop opera¬ 
tors (WLO) along the non-contractible 1-cycle in the y 
direction [T3]. In the thin-torus limit of a general 2-l-lD 
system with TO, the loop order parameters are operators 
that measure the anyon type in the non-contractible loop 
along y direction. The long-range correlation of such or¬ 
der parameter simply reflects the fact that the anyon type 
measured at different x locations should be the same. 
The classical states of this CO are the states with a fixed 
anyon type, which are known as the minimally entan¬ 
gled states (MES) [1]. In our LRMI approache, we can 
choose the LRMI I{D^ ^ describe the TO of the 

2-l-lD system (with 7 the non-contractible 1-cycle along 
the y direction). In the thin-torus limit (with the ratio 
dist(H, B)/Ly effectively infinity), I{D^ x Ty) reduces to 
the LRMI T{D^) that probes the CO in 1-l-lD system 
(see Fig. [T] (a)). Therefore, based on the discussion of 
the LRMI on system with CO, we conclude that, in the 
presence of TO, the LRMI I{D^ x Ty) in the thin-torus 
limit stays finite and display the Shannon entropy form 
for generic states in the Bg and vanishes only on the 
MES states, which were identified above with the classi¬ 
cal states in the quasi-l-l-lD system. 

These results, while initially obtained in the thin- 
torus imit, remain to hold even when Ly/Lx ~ 1. For 
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Figure 1: (a) X{D^ x Ty) are the LRMI between regions 

A and B that wrap around the y direction of the torus. In 
the thin-torus limit, X{D^) reduces to the LRMI X{D^ x Ty) 
that probes the CO in the quasi 1 + lD system, (b) On the 
(edge with the same arrows are identified), the long range 
limit of the pair {Ia,b , Ia' ,b') is constrained by the <S-matrix 



I(D'xTI) 


Figure 2: For the toric model (a), the Ising anyon model(b) 
and the SU{ 2)3 anyon model (c), the permissible regimes of 
{X{D^ X T^),X{D^ X Ty)) are obtained by applying Eq. |^and 
by Monte-Carlo sampling through Bg- 


any non-contractible prime 1 -cycle 7 , we can choose the 
MES basis of the Bg with i labeling the 

anyon types measured by WLOs along the 1-cycle 7 . 
A generic ground state |'0) G Bg can be expanded as 
IV') = Similar to the thin-torus limit, the 

WLO’s following the same type of 1-cycle 7 at different 
locations exhibit non-trivial correlations which is ensured 
by Eq. to reflect itself in the I{D^ x T}^). Using the 
topological field theory and the replica trick, we can show 

m 

X{D^ X T}^) = - ^ log (4) 

which is the Shannon entropy form that agrees with our 
thin-torus analysis and our proposal. Similar to the 
LRMI redefinition of classical states in CO, Eq. [fallows 
us to redefine the MES states as those that minimize the 
LRMI X T.^). As we see here, the LRMI unifies the 
language for the discussion of CO and TO. 

In principle, a gapped 2-l-lD system can possess both 
CO and TO. While the LRMI X{D^) is still a good probe 
of the CO, I{D^ X TJ) receives contribution from both. 
The difference of these two are the contribution from TO. 

Topological uncertainty principle in 2-|-lD 
topological order - We have seen that TO and CO 
are distinguished by LRMI of different regions. There 
is another key feature of TO that is different from CO. 
The loop order parameters in different non-contractible 
cycles generically do not commute with each other. In 
other words, a state with a fixed anyon type in one 
loop is generically a superposition of states with differ¬ 
ent anyon types in the other loop. For the two prime 
1-cycles 7 i and 72 on the torus, their corresponding MES 
basis and {\i )are different and 

related by the transition matrix AI( 7 i, 72 )y = 720 I*) 7 ij 
also known as the modular matrix that carries the defin¬ 
ing data of the TO [51]. For example, if we choose 71,2 
to be the cycles along the x and y direction, Ai(x,y) is 
the modular iS-matrix, which is one of the generator of 
the modular transformation on T^. 

The non-trivial modular matrix AI( 7 i, 72 ) relates the 
LRMFs X T.^\) and I{D^ x T^) through Eq. 

In particular, the 5-matrix constrains the pair {X{D^ x 


Tl),T{D^ X T^) (which is the long-limit of {Ia,b, Ia',b') 
shown in Fig. |R(b)). In Fig. we obtain the permissible 
pair {X{D^ x '7^),I(ZI^ x Ty)) by applying Eq. j^and the 
Monte Carlo sampling through the ground state Hilbert 
space Bg for well-known theories ((a) 2 -l-lD toric code 
model, (b) Ising anyon and (c) the SU{2)3 anyon model). 
We notice that x T^) 7 ^ I{D^ x Ty) generically. 

Furthermore, even though I(ZI^ xT^) and X{D^ xTy) can 
individually take values from 0 to logdim('HG)) the sum 
X T^)+I{D^ ^'^y) always nonzero, reflecting the 
fact that no state in Bg can be the MES with respect to 
both 1-cycles x and y. In fact, by applying the Maassen- 
Uffink inequality |14j. we can obtain 

X T4) +I{D^ X T4) > -21ogmax|M(7i,72),^.| 

'^•’3 

(5) 

This lower bound is indicated by the dashed line in Fig. 

In contrast, the LRMFs in a purely classically ordered 
system do not depend on the topology of the regions and 
thus simultaneously vanish for different region choices for 
the classical states. Therefore Eq [^provides an informa¬ 
tion theoretic measure of the key difference between TO 
and CO. A nonzero lower bound of the sum of LRMI 
in two intersecting pairs of regions measures the non¬ 
trivial commutation relations between the WLOs around 
the two cycles, and guarantees that long-range entangle¬ 
ment is always present in any topological ground state. 
Therefore, we dub Eq. the “topological uncertainty 
principle”. Although the LRMI and topological uncer¬ 
tainty principle appears to be merely a reformulation of 
the well-understood theory of CO and (2-|-l)D TO, it 
provides a new path for understanding TO’s in higher 
dimensions, as we discuss below. 

LRMI in higher dimensional systems with TO 
- Although TO in higher dimensions have been inves¬ 
tigated in some models [I5H37| , the general structure 
of TO in higher dimensions remains an open question. 
The understanding of CO and 2-1-lD TO by LRMI leads 
to a large class of generalizations-the m-membrane con¬ 
densates in which orders in n -I- ID are induced by con¬ 
densing m dimensional membranes (0 < m < n). As 
we have explained above, CO and 2-l-ID TO can be re- 
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Figure 3: In 3+lD, the LRMI of different topologies, includ¬ 
ing (a) ^ X (T^/D^), (b) x , (c) x and 

(d) X {Ty^/D"^), will be considered in the full diagnostic 
of m-membrane condendsate (m = 0,1,2). The non-trivial 
cycles 2 and yz are chosen here as examples. 


spectively viewed as a O-membrane (point) condensate 
and a 1-membrane (string) condensate. Similarly, m- 
membrane condensates with m > 2 have been identified 
in systems with TO in higher dimensions. [nmniiiiiiii 
mUTlEllST] using various of different methods. As we 
shall see, our LRMI approach offers a unified language 
for all TO-membrane condensates in any dimensions. 

For n -f ID systems, we focus on the n-torus T” 
with a linear size L^, and consider two types of LRMI: 

X 7 ^) and x (T™/D'")). Examples 

for 3 -b ID systems are illustrated in Fig. (c) (d). For 
X T™), we consider the LRMI between two re¬ 
gions A and B of the same topology x T™ with 

the T"* part wrapping around a non-contractible prime 
TO-cycle 7 and the linear size r of the D"”™ part fixed 
and satisfying ^ ^ r ^ Lt- The “long range” limit 
of X T™) is defined as dist(A, R) — >■ oo with 

dist(A,R)/LT fixed. I(D"-™ x (T^^/D™)) is defined 
similarly with an extra condition on the (fixed) linear 
size r' of the D™ part: ^ <C r' ^ Lt- With these defini¬ 
tions, we propose 

A gapped n+lD system possesses an m-membrane 
condensate TO (0< m< n) whenI{D'^~"^ xTJ^) > 
y. (j’™//)™)) > 0 for generic states in Hq 
with a given non-contractible prime m-cycle 7 . 

Here we have adapted the formal notation D” x T° = D” 
and Z(D" X iT°/D°)) = 0. 

In the case of n = 2 and m = 1, as we notice that 
X (T^/D^) is topologically equivalent to this pro¬ 
posal is reduced to the proposal for the 2-l-lD TO. In 
3-l-lD, Fig. shows examples of regions with different 
topology relevant to the LRMI: x (T^/D^), 

X X Ty^ and x {Ty^/D'^). Here, the sub¬ 

scripts z and yz represent the non-contractible prime I- 
cycle and 2 -cycle along the 2 direction and the z — y 
plane respectively. As an example, we can show that 
the LRMI indicates the coexistence of 1-membrane and 
2-membrane condensation in the 3-1-ID toric code model 
(or equivalently the Z 2 gauge theory) [13], which agrees 
with the result in Ref. m- In 4-l-lD, two generalized 
toric code models introduced in Ref. [351133 exemplifies 
two fundamentally different 4-l-lD TOs: (1) the coex¬ 
isting 1-membrane and 3-membrane condensates and (2) 


the 2 -membrane condensate m- 

In general, an m-membrane condensate in n -|- ID 
should be equipped with the m dimensional Wilson sur¬ 
face operators (m-WSO) as its order parameter. The 
non-trivial long-range correlation function between m- 
WSOs contributes (by Eq. to all LRMFs x 

T™') (m' > m) and D""™ x (T™"/D™") (m" > m). 
As we notice that the regions x T™ and D"”™ x 

(T™/!?™) contain exactly the same types of WSO (of di¬ 
mension lower than m) but the m-WSO (with dimension 
m). Therefore, the difference between !(£)"“'" x T™) 
and X {Tjp/D'^) only comes from the m-membrane 

condensation. 

Similar to the discussion of 2-l-lD TO, all the TO’s 
signified by non-vanishing x T™) can be viewed 

as the CO in the quasi-(n — m) -b ID system where the di¬ 
mensionality associated to the m-cycles 7 are considered 
as the internal degrees of freedom. Therefore, we gener¬ 
ally expect the Shannon entropy form of x T™) 

which is minimized on the eigenbasis of the WSO’s 
contained in the x T™ region. For two inter¬ 

secting contractible prime cycles 71 ^ 2 , the WSO’s con¬ 
tained in X and x generically 

do not commute, resulting in two distinct sets of eigen¬ 
basis and, hence, the topological uncertainty principle of 
X T™i) -bX(D”-'"2 X r™=). Here mi ,2 denotes 
the dimensionality of the cycles 71 , 2 - In particular, this 
topological uncertainty principle is related to modular 
transformation on T" when mi = m 2 > 0, as it does 
when n = 2 and mi = m 2 = I. 

Discussion - Compared to other entanglement mea¬ 
sures of TO in gapped systems, such as the topological 
entanglement entroDv |401 141] and other similar universal 
terms in entanglement entropy [5311331 - H3| . the LRMI has 
the advantage of being ultraviolet finite. All subtleties 
and ambiguities in entropy calculation due to entangle¬ 
ment of short-scale degrees of freedom [331 135337] cancels 
in LRMI. As the universal terms of the entanglement en¬ 
tropy also appears in gapless systems [331SS] 135353] , it is 
interesting to also generalize the LRMI proposal to these 
cases. In particular, different from disks, the size of the 
non-contractible regions have to grow with the system 
size, resulting possible extra contributions to the LRMI 
from the gapless mode. The criteria of ordering needs to 
be modified accordingly. 

Ref. |5I| proposed the non-conventional orders in¬ 
duced by the spontaneously broken m-form global sym¬ 
metries in a general dimension. We believe that the 
LRMI not only captures such non-conventional orders as 
m-membrane condensations but also provides more infor¬ 
mation because the long-range correlation does not need 
to be associated with a group structure. As is exemplified 
by 2-blD TO, SU{2)k TO states share the same 1-form 
symmetry [5T] but have distinct LRMI and topological 
uncertainty relations. It would be important to study the 
detailed connection between the proposed m-membrane 
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condensate and the m-form global symmetries. 
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Supplementary material 

This supplementary material contains 4 sections. In the first section, we review the definition of the 2+lD toric 
code model and study its thin-torus limit, demonstrating the connection between the 2+lD topological order and the 
conventional order in 1+lD. The second section provides a derivation of the long-range mutual information (LRMI) 
Eq. 4 for the 2-torus and generalize this equation to a general 2-manifold. In the third section, we study the LRMI 
in the 3-l-lD toric code model. We show that this model exhibits both 1-membrane and 2-membrane condensations. 
In last section, we study the LRMI in 2 different 4-l-lD toric code models with one of them hosting a coexisting 
1-membrane and 3-membrane condensations while the other exhibits the condensation of 2-membranes. 


THIN-TORUS LIMIT OF THE 2-|-lD TORIC CODE MODEL 


In this section, we will show how the Z 2 topological order (TO) of the 2-l-lD toric code model reduces to the 
conventional order (CO) of the quasi 1-t-lD system in the thin-torus limit. The 2-l-lD toric code model is defined on 
a 2D square lattice with the degrees of freedom on the links and the Hamiltonian given by 


Htc = “ X! “ X! 


(SI) 


where v and p labels the vertices and the plaquettes respectively (see Fig. 
plaquette term Bp are given by 


SI (a)). The vertex term Ay and the 


Ay = Y[af, Bp = Y[c 


(S2) 


l^V 


iGp 


where the index I labels the links that connect to the vertex v or belong to the plaquette p and denote the Pauli 
matrices. All the vertex terms Ay and plaquette terms Bp commute with each other, making this model exactly 
solvable. A ground state \tj}) of H^c should satisfy 


Ay\tjj) = Bpl'lp) = |' 0 ), 


(S3) 


for any vertex v and plaquette p. The ground state of L7 tc on an S'^ is non-degenerate. However, there are 4 
degenerate ground states on the torus This 4-fold ground state degeneracy is a signature of the TO. The 4 
degenerate ground state can not be distinguished by local operators, while the non-contractible Wilson loop operators 
(WLO) act non-trivially among them. The WLO are given by 


IPe(C) = n' 

lec 


Wr 


AC )=n 

lec 


(S4) 


where C and C are both non-contractible closed path on the T^, one along the links of the lattice and the other along 


the links of the dual lattice (see Fig. SI (a)). It is easy to show that 


[We(C), Htc] = 0, [W„(C'), Htc] = 0. 


(S5) 


When acting on the ground state Hilbert space Hg, these WLO’s We{C) and Wm(C) do not depend on the detail of 
the path C and C' , but rather only on the type of non-contractible 1-cycles of the path C and C traverse along. 
In Fig. SI (a), we have chosen the path C and C to follow the non-contractible 1-cycle along the y direction. In 


this case, they commute with each other and each of them have eigenvalues ±1. The 4 simultaneous eigen states of 
We{C) and WmiC) with 4 different combination of the eigenvalues serve as a basis of Hg on the torus These 
eigenvalues in fact label the different anyon types in the TO threading through the non-contractible 1-cycle along the 
y direction. Therefore, as is proposed in Ref. [T] , this basis coincide with the basis of the minimally entangled states 
(MES). 

Now, we consider the thin-torus limit of the 2 -|- ID toric code model. Because the 2 -|- ID toric code Hamiltonian 
consists of only commuting terms, the correlation length of the system is in fact 0. Therefore, we can, without loss of 
the topological nature of the system, take the thin-torus limit with the circumference in the y direction Ly = 1 in the 


lattice unit (see Fig. SI (b)). In this limit, the model is reduced to two decoupled Ising chains with the Hamiltonian 


77idTC = — ^ 


<^1+1 


2 G Z 


i^Z 




(S6) 







Figure SI: (a) The 2+lD toric code model is defined on a 2D square lattice with the degrees of freedom (black dots) on the 
links. In its Hamiltonian Htc, the vertex term acts on the 4 links (red) connected to the vertex v, while the plaquette term 
Bp acts on the 4 links (blue) on the boundary of the plaquette p. The Wilson loop operator We{C) and Wm{C') are defined as 
a product of spin operators along the path C (orange) and C' (green), (b) In the thin-torus limit where Ly = 1 in the lattice 
unit, the 2+lD toric code model is reduced to a spin chain shown in the figure. The WLO (orange and green), 114(C) and 
Wm{C'), reduce to single spin operators and the local order parameter of the spin chain. 


The subscripts i and f + 1/2 with i G Z label the links along the x and y direction respectively. Obviously, the system 
has a global Z 2 x Z 2 symmetry and the ground states possess the conventional order as 1+lD system. In this limit, 
the Wilson loop operators We{C) and Wm{C') are reduced to single-spin operators trf and cr /^^^27 which are exactly 
the order parameters of the 1-l-lD system. Therefore, the MES of the 2-l-lD TO are identified with the classical states 
in the 1-l-lD conventional order. 


LONG-RANGE MUTUAL INFORMATION FOR 2-|-lD TOPOLOGICALLY ORDERED SYSTEMS ON 
THE 2-TORUS AND OTHER GENERAL 2-MANIFOLDS 


In the main text, we present Eq. 4 for the long-range mutual information (LRMI) in 2-l-lD topologically ordered 
system on a 2-torus . In this section, we provide a field theoretic derivation of Eq. 4 and generalize it to a general 
2-manifold. First, we review the calculation of the entanglement entropy for 2-l-lD topological states in a general 
setting by using the replica trick [5]. Consider a general 2-l-lD topological system with the degrees of freedom given 
by the field (p{x) and the (Euclidean) field theory action given by Se[‘^(x)]. In the limit where the physical correlation 
length ^ —7> 0, the (unnormalized) ground state wave functional fl'[:^(x)] on a closed 2-manifold Af can be written as 
a path integral on any 3-manifold B such that dB = M with Lp{x) as the boundary condition: 

^[(^(x)] = f D[^'(y)]e-®^[^'(^)1. (S7) 

The (unnormalized) reduced density matrix of a region X on the 2-manifold M can also be expressed in the path 
integral form. Consider two copies of B denoted as B^ and B^ . The corresponding regions on their boundaries are 
denoted as X^ and X"^. We partially glue together the 3-manifolds B^ and B'^ along dB^ jX^ and dB"^ jX'^^ creating a 
new 3-manifold 77. with its boundary dK, = X^ U X^. The (unnormalized) reduced density matrix can be written as 


pxMx),(p'{x)]=Tr^/x\'^){'i’\= f D[(p"(y)]e 

Jy^n, ¥3"(i/)|^i=v(i/), v"(y)\^ 2 =^'iy) 

The entanglement entropy of the region X can be obtained as 


Sx = - lim — 
n-s-i dn 


( Trp" \ 


(S8) 


(S9) 


We can express Trp)(^ as a partition function Z(77„) on the replica manifold 77„. 77„ is obtained from gluing n copies 
of the 3-manifold 77 where the X^ part of the boundary 977 of the fth copy is identified with the X^ part of the 
boundary 977 of the i -I- 1th copy (with i = l,2,...,n and the identification n -|- 1 ~ 1). Therefore, we obtain the 
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entanglement entropy as 


Sx = — lim 


1 dn 


( ZjT^n) \ 
VZ(7^l)"/ ■ 


(SIO) 


Derivation of the LRMI Eq. 4 on 2-torus 


Eq. 4 applies to the LRMI I^{D^ x T^) for any non-contractible prime 1-cycle 7 on T^. Since any non-contractible 
prime 1-cycle 7 can be mapped to each other by the modular transformations of the T^, it is sufficient to derive Eq. 
4 just for the LRMI ^ T.^) for a specific 1-cycle 7 which is chosen here to be the non-contractible prime 1-cycle 

along the y direction. For this purpose, we need to calculate the entanglement entropy Sa, Sb and Saub for regions 


A and B shown in Fig. S2 (d). Ref. [2] has already presented the derivation of Sa and Sb which we will review for 
the sake of completeness. Then, we generalize this calculation to Saub and derive the LRMI I^{D^ ^ )■ 

Now we restrict our discussion to the ground states on the 2-torus N = with the corresponding 3-manifold 
B — X S^ a solid torus. We will assume that a 2-l-lD topological system always admit a Chern-Simons theory 
description. This assumption allows us to construct different degenerate ground state on the by inserting different 


WLO’s into the path integral on ,B[3] (see Fig. S2 (a)): 








(Sll) 


This wave functional exactly represents the MES state \i)y (defined in the main text) with ith type anyon 

threading through the non-contractible 1 -cycle along the y direction. 

we will calculate the entanglement entropy Sa of the x region A on 


S2 


(a)). 


Following Ref. [2], we will calculate the entanglement entropy Sa of the x Ty region A on (see Fig. 

We first focus on the entanglement entropy Si^A on the MES state \i)y (or equivalently 4'i[(/3(a;)]). The solidMorus 
B = D'^ X S^ from which 4'i[(/?(a;)] is obtained is topologically equivalent to 2 3-disks connected by two “bridges” 
S^ X (see Fig. S2 (b)). We denote the reduced density matrix of \i)y on the region A is denoted as pi^A- The 
quantity Tv{pi A)^ in Eq. [^is equivalent to the partition function on the replica manifold Tin obtained from gluing 
2n copies of the solid torus B = D'^ x S^ together. The gluing process is done by identifying the boundary of these 
2n solid tori according to the dotted blue lines shown in Fig. S2 (c). In Fig. S2 (c), the 3-disks on the upper row 


are glued together into an 3-sphere S^ with 2n 3-disks punctures. The same applies to the 3-disks on the lower 
row. The 4n “bridges” are pairwise glued into 2n “tubes” of the topology S^ x D^, connecting the two 3-spheres 
S^ through their punctures. Therefore, the replica manifold Tin is a manifold with 2 3-spheres S^ connected by 
2n 
the 


“tubes”. The partition function on Tin can be evaluated through surgery Hig. As is shown in Fig. |S3[ each of 
“tubes” that connects the two 3-spheres S^ can be disconnected into two pieces with the two open ends capped 
off by two 3-disks and the Wilson loops reconnected accordingly. The partition function of the manifolds before 
and after this surgery process differ by a factor of (iSoi)”^, where S is the topological 5-matrix, i is the anyon type 
the WLO Wi carries and 0 represents the trivial anyon. After disconnecting all the “tubes” in Tin, we obtain two 
disconnected S^ each carrying a WLO W^, the partition function of which is given by |5oip. Therefore, we have 


Tr(p"^) = 


(S12) 


Applying Eq. |S9[ we obtain 


= =21og|5o,|, (S13) 

n^l dn \ ( LTpi^A) / 

which is the topological part of the entanglement entropy HIS]- The usual “area law” term in the entanglement 
entropy vanishes in this calculation because the application of the Chern-Simons theory effectively projects out the 
higher energy physics, leading to a ^ = 0 correlation length. 

Now we can caculate the entanglement entropy S^^a for a generic state |'0) in the ground state Hilbert space TLg on 
T^. We can expand the state \tp) in the MES basis: \'ij)) = '0y,i|*)y Es reduced density matrix p^^A on the region 
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A then takes a block diagonal form: 






. Therefore, the entanglement entropy is given by 


S,I,,A = ^ ^ log IV'y.iP 

i i 

= XI 2|V'y.z P log |5o*I - X \'^vA^ log \Av,^ 


(S14) 


With the same calculation, we can show that the entanglement entropy for the region B is identical to S^^a- 
In the following, we will generalize the calculation of studied in Ref. [2] to the entanglement entropy S^^aub 
for the regions AU B shown in Fig. [^(d). Again, we first focus on the MES state \i)y ( or 'I'i[(/j(x)]). The solid 
torus B = X from which the wave functi ona l d>i[(p(a:)] is obtained is topologically equivalent to 4 3-disks 
connected by four “bridges” x (see Fig. S2 (e)). The reduced density matrix of \i)y on the region A U R is 


denoted as pi^AuB- Ti'(p"yiuB) 1 ^® identified as the partition function on the replica manifold obtained from 
gluing 2n copies of the solid torus B = x together in the way indicated in Fig. S2 (f). Similar to the discussion 
on Tin, we can view the replica manifold Tin as 4 3-spheres S^, each with 2n 3-disks punctures, connected in 
a cyclic order by in total 4n “tubes” x D^. Through surgery, we can obtain the partition function on Tin and 
therefore Tr(p”^^jg) as 




Tthe entanglement entropy is then given by 


Si,AUB = 4 log |5oi|. 


(S15) 


(S16) 


For a generic state \ip) = '^i''Py,i\i)y fh® ground state Hilbert space, the reduced density matrix pp^AuB on the 
region AU similar to also takes a block diagonal form. Therefore, the entanglement entropy is given by 


^ ^ Si^AUB ^ ^ l*^g 

i i 

= X'^l^y^l^lOgl^Oil - X \Av,if^Og\^y,i\^. 


(S17) 


Since we have taken the limit of vanishing correlation length, the LRMI x Tj) is directly given by the mutual 

information between regions A and B : 


X Ty ) — S^^A 4“ Bip^AUB — ^ ^ fog \'ipy,i 


(S18) 


which is exactly the Eq. 4 of the main text. In the presence of a finite correlation length the mutual information 
between A and B has an extra contribution from the local correlations with expected form ^ In the 

definition of LRMI, we consider the limit with dist(A, B)/Lx,y finite and dist(A, B) -A oo. The contribution from the 
loca l cor relations vanishes in this limit. Therefore, the result of x Ty) is still valid. It is worth mentioning that 

Eq. S18 allow us alternatively define the MES basis {|*)j/} as the states that minimizes the LRMI I^{D^ x Ty). 


Derivation of the LRMI X,p{D^ x Tj) on a general 2-manifold 


In fact, the derivation of the LRMI I^{D^ x Tj) in the previous subsection can be generalized to a general 2- 
manifolds Af with any non-contractible prime 1-cycles 7 . We can always deform the 2- mani folds JV such that the 
neighborhood containing the two disconnected x Tj region A and B coincides with Fig. 
the rest of the manifold Af outside of this neighborhood. Let us first consider a state \i,u., 
i threading through the 1-cycle 7 . The wave function on U is specified by the label u, 


S 8 


(a). Here 7/ represents 


with a fixed anyon type 
The path integral on B that 

generates this state \i,Ui) should contain a WLO as is shown in Fig. S4 (a). The entanglement entropy Si^a of 
the state \i,Ui) on the region A can be calculated through Tr(pi^yi)". Tr(^^) which is equivalent to the partition 
function on the replica manifold Tin obtained from gluing 2n copies of B together in the way indicated by Fig. S4 
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Figure S2: (a) The MES \i)y on the with the ith type threading threding through the non-contractible 1-cycle along the 

y direction can be generated by the path integral on the solid torus B = x with the WLO inserted. The x Ty 
region A of the is depicted, (b) shows a 3-manifold that is topologically equivalent to the solid torus B = x . (c) 
shows the identihcation (indicated by the blue dotted lines) in the process of gluing 2n copies of B into the replica manifold 
TZn- (d) depicts the two separated regions A and B on the T^, each with the topology of x Ty. (e) shows a 3 manifold 
that is topologically equivalent to the one shown in (d). (f) shows the identification (indicated by the blue dotted lines) in the 
process of gluing 2n copies of B into the replica manifold TZn- 




Figure S3: The surgery procedure relates, by a factor of the partitions on two different closed 3-manifolds. The 3-manifold 
on the left has R x 5^ structure (a “tube” structure) in a certain neighborhood and the rest of the 3-manifold is denoted as 
KL. The 3-manifold on the right is obtained from the previous one by switching the “tube” structure to ((R~ x S'^)/({0“} x 
S^)) U ((R”*" X •S'^)/({0’'"} X S^)) and reconnecting the WLO’s without changing KL. This procedure effectively disconnects the 
“tube” structure into two pieces with the two open ends capped off by two 3-disks . 



(b). Notice that, by performing the surgery procedure along the green dotted line in Fig. S4 (b), we can isolate the 
contribution from the 2n copies of the lA part (including its interior) of the replica manifold TZn and obtain 


(S19) 


where Jyy(|f,rti)) is the contribution from one copy of the lA (and its interior), which depends on the quantum state 
\i,Ui). Regardless of the specific form of .Fu, we always have 


MpIa) 

Tr(p,.A)" 


and, thus, the entanglement entropy is given by 


|5o*| 


2(l-n) 


= 2log \Soi\, 


(S20) 


(S21) 


which depends only on the anyon type i but not the other wave functions label Ui- We 
is exactly the same as result Eq. S13 on T^. The same result applies to the suregion 
derivation can be done to obtain the entanglement entropy Si^AuB of the state \i,Ui) 


also recognize that this formula 
B: Si^B = 21og|iSoi|. A similar 
on the region A U R: 


Si,AvjB = 4 log |5oi|, 


(S22) 


which again coincides with the result on the and only depends on the anyon type i. 
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(a) 




Figure S4: (a) shows the manifold M with a neighborhood containing the region A and B and with the rest of manifold 

denoted as U. (b) shows the identification (indicated by the blue dotted lines) in the process of gluing 2n copies of B into the 
replica manifold TZn- The surgery procedure along the green dotted lines allows us to isolate the contribution from U and its 
interior to the partition function. 


A generic state |^/>) on Af can always be expanded in the basis with fixed anyon type threading through the non- 
contractible 1-cycle 7 : \'il)) = . Similar to the case oi N = T^, all three reduced density matrix 

and P' 4 ,,aub are block diagonal (with each block labeled by the anyon type), following which we have 

X T^) = + S,,B - S^^aub) - log 

i i 

= “Xl (S23) 

i 

This result again follows the Shannon entropy (SE) form on the probability distribution {|' 07 ,ip}i of the quantum state 
IV') on the anyon type threading through the non-contractible 1-cycle 7 . The LRMI x Ty) is still minimized 

on the MES state on a general 2-manifold Af 


LONG-RANGE MUTUAL INFORMATION IN THE 3-|-lD TORIG CODE MODEL 

In this section, we will study the long-range mutual information of the 3-l-lD toric code model on the 3-torus T^. 
The 3-l-lD toric code model is defined on a cubic lattice with the degrees of freedom on the links and the Hamiltonian 
given by 


^^3DTC — — ^ ^ Bp, 


(S24) 


where v and p labels the vertices and the plaquettes. The vertex term Ay and the plaquette term Bp in the Hamiltonian 
^^ 3 DTC are given by (see in Fig. S51 


Ay=Y[ar, Bp = Y[, 


7 I 


(S25) 


iGv 


iGp 


where the index I labels the links that connect to the vertex v or belong to the plaquette p. All terms in the Hamiltonian 
commute with each other rendering this model exactly solvable. Any ground state \tjj) of the Hamiltonian should 
satisfy 


Ay\ij)=Bp\^) = \^), (S26) 

for any vertex v and any plaquette p. The ground state is non-degenerate on a 3-sphere. However, the ground state 
degeneracy on 3-torus is 8 -fold indicating the existence of TO. Within the ground state Hilbert space 'Hq of the 3 
torus T^, the WLO Wx,y,z and 2-WSO’s Wxy,yz,zx acts non-trivially. (Here we’ve adopted the notation “2-WSO” 
from the main text.) Wx is a product of along the links which form a path that winds around the non-contractible 
prime 1 -cycle along the x direction. Wy^z are defined similarly. Wxy is a product of on the links whose dual 
plaquettes form a 2D surface wrapping around the non-contractible prime 2-cycle along the x — y plane. Wyz,zx are 
defined similarly. 

In the following, we will follow the proposal in the main text and show that the 3-l-lD toric model exhibits 1- 
membrane and 2-membrane condensation by computing the following LRMI on the 3-torus T^: X{D^), I{D^ x T^), 
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Figure S5: The vertex term Av and the plquette terms Bp are depiced in this figure. 


X{D^ X Ty^) and X{D^ x (Ty^/D^)). These calculations are done for a generic ground state \tp) in the ground state 
Hilbert space Hg on the T^. It is worth mentioning that Ref. [TB] also pointed out that the 3+lD toric model can 
be viewed as a condensate of “strings” and “membrane” (which correspond to 1-membrane and 2-membrane in our 
proposal). However, our definition of condensate is completely different from, and should be more general than the 
study presented in Ref. [B]. Before getting into the detailed calculations, we first set up the convention for the choice 
of regions and terminology. As a convention, we define a region to be a collection of cubic cells such that if a cubic 
cell belongs to this region, so do all of its vertices, links and plaquettes. Also, if a vertex v, a link I or a plaquette 
p belongs to the region, there must exist a cube in this region such that this cube contains v, I or p respectively. A 
vertex is defined as a boundary vertex of a region if the vertex v belongs to the region, but there exist links outside 
of this regions that are connected to v. 


Calculation of the LRMI I{D^) 

We consider two separated regions A and B. The entanglement entropy and S^^agb associated to 

the regions A, B and A U R can be obtained by the Schmidt decomposition of First, we consider the region A. 
For each boundary vertex v G dA, we define the boundary vertex operator A„ = ^ a product of cr“’s on 

the links that both belong to the region A and are connected to the vertex v. A boundary condition on dA is defined 
by the set which is a choice of eigenvalues ± 1 for all the boundary vertex operators Ay. The ground 

state IV') admits a Schmidt decomposition: 

IV') (S27) 

i 

with the Schmidt state j) defined on the region A satisfying the boundary conditions 

Vu G dA, (S28) 

and the “ground state conditions” 

Vu G {A/dA) 

TJpIV'fa.}) = IV'(1})> VpgA. (S29) 

Since the number of equations above and the number of degrees of freedom contained in the region A, the Schmidt 
state j) is uniquely determined by these equations, which confirms that the one-to-one correspondence between 

the Schmidt states and the choices of boundary conditions. Notice the identity Wv^qa^v = Y\v^(A/dA) A 
boundary condition is permissible only when 

a, = 1. (S30) 

v^dA 

The number of permissible boundary conditions is then given by with \dA\ the number of boundary vertices. 

Also notice that, for the plaquettes p that share links with but don’t belong to the region A, the action of Bp leaves 
the ground state l"^) invariant, but maps the Schmidt state with one boundary condition to another. The Schmidt 
coefficients of boundary conditions that are related by these plaquette operators Bp therefore have to be the 

same. Since dA doesn’t contain any non-contractible 2-cycle of the T^, all the boundary conditions on dA are related 
to each other, rendering the Schmidt coefficients A{a„} a constant. We can then obtain the entanglement entropy of 
the region A: 


R^,A = (|aA|-l)log2. 


(S31) 
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The same derivation applies to the region B and results in = i\dB\ — 1) log 2. For the region AU B, its two 
disconnected part A and B leads to two independent constraints on the permissible boundary conditions. The number 
of permissible boundary conditions on d{A\J B) is then given by Again, since the Schmidt coefhcients 

take a constant value, the entanglement entropy oi A\J B is given by. 

5v,,aub = (|5(AuB)|- 2)log2. (S32) 

Notice that the |9(AU B) \ = \dA\ + \dB\. Since the 3+lD toric code model has 0 correlation length, we can directly 
obtain the LRMI I{D^) as 

I{D^) = = 0, (S33) 

which indicates the absence of CO in the 3+lD toric code model. 


Calculation of the LRMI I{D^ 'xT}) 


We consider two separated x Tl regions A and B. The method of Schmidt decomposition of the ground state 
\'(p) is still valid. We first consider the region A. Due to the existence of a non-contractible 1-cycle, namely the 
non-contractible 1-cycle along the z direction, in the region A, the Schmidt states are only labeled by the permissible 
boundary conditions {a„}. They are also eigenstates of the WLO Wz and, therefore, carries the labels Wz = ±1 which 
denotes their corresponding eigenvalues Wz = ±1- Therefore, we can write the Schmidt decomposition as 

\'^) (S34) 

i 

with the Schmidt states satisfying 

Vu e (A/dA) 

Vp G A 

I V'fa„ },»,)■ (S35) 

One can show, by matching the numbers of the equation above with the number of degrees of freedom contained 
in the region A, these equation uniquely determine the Schmidt state \4’fa } m )■ Since the boundary dA does not 
contain any non-contractible 2-cycles of the T^, all the boundary conditions are related to each other, rendering the 
Schmidt coefficients only a function of their Wz label, i.e. A{a„},t«^ = The values of is fixed (up to a U{1) 

phase) by the expectation value of the WLO Wz on the ground state \tp): {Wz)ip = • Hence, we can obtain 

the entanglement entropy S^^a as 


VA = (|9A|-l)log2 


2 ^2 


1 + {Wz)^ 1 + {Wz)^ 

2 2 


(S36) 


The same derivation for the region B shows that = {\dB\ — 1) log2— log 

For the region AUH, the Schmidt state are still labeled by their boundary conditions on d{A US) and the eigenvalue 
Wz of the WLO Wz- Similar to the case of X{D^), the number of permissible boundary conditions is again given by 
2 |a(AuB )|-2 entanglement entropy S,p^AuB is then given by 


O OM o l-{Wz)^ ^+{Wz)^, l+{Wz)^ 

S^,AUB = (|9(AUH)| -2)log2---^ log---^---^ log--- 


and the LRMI X{D‘^ x T^) is given by 

X{D^ X T^) = S^^A + — «S'y,_AUB = — 


1 - (Wz)^ , 1 - {Wz)i, 1 + {Wz)i, , 1 + (Wz)^ 


log 


log 


(S37) 


(S38) 


which is the SE of the classical probability distribution of the state \'tjj) on the two eigenvalues Wz = ±1 of the WLO 
Wz- As is proposed in the main text, this non-trivial value of X{D^ x Tl) signifies the 1-membrane condensation in 
the 3-l-lD toric code model. 
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Calculation of the LRMI T{D^ x (Ty^lD’^)) 


We consider the two separated x {Ty^/D'^) regions A and B. Each of the regions A and B contains two non- 
contractible 1-cycles in the y and z directions along which the WLO’s Wy and can be constructed. Similar to the 
previous case, the Schmidt states rely not only on the boundary conditions, but also their eigenvalues Wy,Wz = ± 
with respect to the WLO’s Wy and Wz. Due to the non-existence of non-contractible 2-cycles of contained in dA, 
dB or d{A U B), the Schmidt coefficients do not depend on the boundary conditions, but only their Wy^z labels. A 
calculation similar to the previous case leads the the conclusion that 


%^ = (|aA|-l)log2- ^ 

'Ye{y,z} 

S^^B = {\dB\-l)\og2- ^ 

76{y,z} 


2 ^2 2 ^ 2 ) 


1 - log 1 - I l + + 


5^.^uB = (|5(AuB)|-2)log2- ^ 

76{!/.z} 


I-{w,)^ I i + + 

2 ^ 2 


2 2 


r 


(S39) 


The LRMI I{D^ x (T^ /D"^)) is then given by 


I{D^ X iTlylD^)) = - S^^aub = - 

76{y,z} 


2 ^22 ^ 2 J 


(S40) 


which coincides with the SE of the classical probability distribution of the state \ip) on the 4 combination of eigenvalues 
{wx,Wy) of the two WLO’s Wy^z- 


Calculation of the LRMI I{D^ x T; 


We consider the two separated x Ty^ regions A and B. For the region A, the Schmidt decomposition of the 


ground state tp, as we discussed in the previous case, is still given by 

^ ' ^{a^},Wy,W^ \'4’{ay},Wy,W;^)\'^{dy 


(S41) 


where {d„} denotes the boundary conditions and Wy^z = ±1 correspond to the eigenvalues of the WLO’s Wy^z- In 
contrast to previous cases, due to the existence of non-contractible 2-cycles which is the one wrapping around the 
y — z plane, all boundary conditions are not completely related to each other. To be more precise, the boundary dA 
of the region A consists of two disconnected parts (9A)i^2, each with the topology of Ty^. The action of Bp terms of 
any plaquette p does not change the value of IlvelOA)! (or equivalently nt 2 e(aA )2 I^ II I® easy to show 

that n? 2 e(aA)i = Wyz- Therefore, we can identify the n 22 e(aA)i with the eigenvalue Wyz of the 2-WSO Wyz, 
which can take values ±1. Wyz provides a label for the topological sectors of the boundary conditions within which 
all the boundary conditions are related to each other. The Schmidt coefficients X{dy},w^,wy can then be rewritten as 


'{5,1, } 


— An 


-n 


^G(dA)i 


(S42) 


which are then determined (up to a C/(l) phase) by the expectation values {Wyz)if,, {Wy)^ and (IFz)^. With the 
values Wyz, Wy and Wz fixed, the number of permissible boundary conditions are given by Therefore, the 

entanglement entropy is 


S'V’.A = (|9A| - 2) log2 


E 

'Ye{vz,y,z} 


2 ^2 


1 + {W,),P „1 + (IF2)2/,^ 

2 2 


J- 


(S43) 
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Similarly the entanglement entropy and S^^aub are given by 

'iniEiklog 1 ^ 1 + {W^)^ 1 + {W,)^ \ 


S^,B = i\dB\-2)\og2- 

1&{vz,Va} 


S„„. = (|a(AUB)|-4)I„g2- i: (i^,„gl^ + l±^,„gl±|!3i), (S44) 


1&{vz,VA} 


Now, we obtain the LRMI X{D^ x Ty^) is then given by 

1(0^ X Tly) = S^^A + S^,B - S^^aub = - Y. 

7e{yz,y,z} 


1 - (m)^ , i-(iRg)^ 1 + {W,)^ 1+JMYt] 

2 2 2 2 


I ? 

(S45) 


which coincide with the SE of the classical probability distribution of the state {ip) on the 8 combination of eigenvalues 
{wyz, Wy, Wz) of the 2-WSO Wyz and the two WLO’s Wy^z- Generically, the LMRI T{D^ x Ty^) differs from the LRMI 
I{D^ X [Ty^/D"^)). Therefore, we conclude that the 3+lD toric code model exhibits the 2-membrane condensation. 


LONG-RANGE MUTUAL INFORMATION IN THE 4-|-lD TORIC CODE MODELS 

In this section, we will discuss the behavior of the LRMI in two different types of 4-l-lD toric code models [71 [S]. Using 
the LRMI diagnostics, we will show that one of them hosts coexisting 1-membrane and 3-membrane condensation and 
the other exhibits the 2-membrane condensation. To be more specific, we will consider the a generic ground state \ip) 
in the ground state Hilbert space 'Hq on the 4-torus T'^. The relevant LRMI here are x T^), I{D^ x T^y), 

Z(U2 X (r4/D2)), Z(U1 X T4J andI(Di X (r4,/D3)). 

Let’s first set up some conventions and terminology. The 4 spatial dimensions are denoted as x, y, z and u. Both 
models are defined on a 4D hypercubic lattice. We will then refer to a vertex as a 0-cell, a link as a 1-cell, a 2D 
plaquette as a 2-cell, a 3D cube as a 3-cell and a 4D hypercube as a 4-cell. If an m-cell Cm contains an n-cell c„ 
(n < m), we interchangeably denote this relation as Cm € c„ or c„ S Cm- For our LRMI diagnostic, we definite a 
region to be a collection of 4-cells such that if a 4-cell belongs to this region, so do all of its vertices,links, plaquettes 
and 3-cells. Also, if a 0-cell u, a 1-cell I, a 2-cell p or a 3-cell c belongs to this region, there must exist a 4-cells in the 
region such that this 4-cell contains v, I, p or c correspondingly. A vertex v is defined to be on the boundary dA of a 
region A if v G A and there exists a link I such that the link connects to the vertex v, namely I G v, but I ^ A. A link 
I is defined to be on the boundary dA of a region A if I G A and there exists a plaquette p such that the plaquette 
contains the link I, namely p G I, but p ^ A. 


4-|-lD toric code model with 1-membrane and 3-membrane condensation 

The 4-l-lD toric model that exhibits the 1-membrane and 3-membrane condensation is defined on the 4D hypercubic 
lattice with the degrees of freedom on the 1-cells (links) with the Hamiltonian given by 

= - E E (S 46 ) 

0-cell V 2-cell p 

where the vertex operators A„ and the plaquette operators Bp are given by 

A,= n af, Bp= n (S47) 

1-cell iGv 1-cell iGp 


All the terms in commute with each other, rendering this model exactly solvable. The ground state of 

non-degenerate. However, there is 2‘^-fold ground state degeneracy on the 4-torus T^, which is 
a signature of topological order. Within the ground state Hilbert space Ha on the WLO’s, including Wx, Wy, Wz 
and Wu, and the 3-WSO’s, including Wxyz, Wyzu, Wzux and Wuxy, act non-trivially. The subscripts of the WLO’s 
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and the 3-WSO’s represent the non-contractible prime 1-cycles or 3-cycles of the these operators are constructed 
along. Each of them has eigenvalues ±1. All of the non-trivial commutation relations among them are given by 

W^Wy.u = -Wy,uW^, WyW.uX = -W,uxWy , W.W^^y = - IE, , UUd W^y, = -W^y.W^- (S48) 

Now we present a general discussion for the region of T'^ with all possible topologies relevant to the LRMI di¬ 
agnostics. On a surbregion A, we can define, for every boundary vertex v € dA, a boundary vertex operator 
Ay = Jli.ceii I lev & leA^i- ^ choice of eigenvalues d„ = ±1 for each of the boundary vertex operator Ay will be 
referred to as a boundary condition {dy}y^gA on dA. For all the relevant cases, the following Schmidt decomposition 
applies: 

1^) “ (S49) 

i 

where is a set of choices of the eigenvalues w.y^ = ±1 for all the WLO’s contained in the region A. 

The boundary conditions are all related to each other unless non-trivial topological sectors exists. For every non- 
contractible 3-cycles 73 of the contained in the region A, there is a label = ± for the topological sectors of the 
boundary conditions. This label w.y^ correspond to a choice of the eigenvalue of the associated 3-WSO Only 

the boundary conditions with the same topological sector are related to each other. Therefore, the values of Schmidt 
coefficients only depends on the labels {^ 73 } and the When the non-contractible 1-cycles or 3-cycles does not 

exist, the corresponding labels drop out automatically. Also, for all the relevant cases, we notice that the numbers of 
permissible boundary conditions satisfy that if the region A is a union of two disjoint regions (say A^ and A^), the 
number of permissible boundary conditions on dA is a product of the number of permissible boundary conditions on 
9A^ and that on dA^. With these general understanding of the Schmidt coefficients, we can derive the all the relevant 
LRMI. The result are summarized as follows: 


= 0 , 

I{D^ X T^) = - ^ 






E (l^logl^Tl±f^logl±f^), 


76{a:.y} 

I{D^ X = - E 

'ye{x,y,z,xyz} 

I{D^x{TlyjDA = - E 

~ie{x,v,z} 




log ■ 


J 


2 ^2 2 ^ 2 j 


(S50) 


The generically non-vanishing values ofX(D^ xT^j)—andZ(i4^ xT^y^)—I{D^ x {T^y^/DA signify the existence 
of 1-membrane and 3-membrane condensations. In contrast, the fact that 0 = I{D‘^) = X{D‘^ xT‘^y)—I{D^ x {T^y/D'^)) 
indicates the absence of the 0 -membrane and 2 -membrane condensations. 


4+10 toric code model with 2-membrane condensation 

The 4-l-lD toric model that exhibits the 2-membrane and 3-membrane condensation is defined on the 4D hypercubic 
lattice with the degrees of freedom on the 2-cells (plaquettes) with the Hamiltonian given by 

E E (S51) 

1-cell I 3-cell c 

where the link operators A( and the cube operators R' are given by 

A= n n (S52) 

2-cell pGl 2-cell pGc 
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This model is exactly sovable, since All the terms in commute with each other. The ground state of 

is non-degenerate. However, there is 2®-fold ground state degeneracy on the 4-torus which is a 
signature of topological order. Within the ground state Hilbert space 'Hq on the T^, the 2-WSO’s, IT^, 

1T“„, Wy^ and {a S {e,m}), act non-trivially. The subscripts of the 2-WSO’s represents the non-contractible 
prime 2-cycles they are constructed along. The superscripts e and m represent the electric WSO’s and the magnetic 
WSO’s respectively. Each of them has eigenvalues ±1. All of the non-trivial commutation relations among them are 
given by 


W® = 

* * xy * zu 

w® W’" = 

'* XU yz 


* * zu* * xy ? * * yz * * XU 

**yz**xui * * yu * * zx 


* * XU * * yz^ * * zx * * yu 

-wzwL. w:,wz = 


* * yu * * zx"! 

* *xy* *zu 


(S53) 


Similar to the previous case, we will present a general discussion for the region of with all possible topologies 
relevant to the LRMI diagnostics. On a region A, for every boundary link I G dA, we can define the boundary link 
operator AJ = J| 2 -ceii p, pel & peA^p- ^ choice of eigenvalues d[ = ±1 for each of the boundary link operator AJ 
will be referred to as a boundary condition {dj};gg^ on dA for this model. For all the relevant cases, the following 
Schmidt decomposition applies: 




A 

afl.f™? 


(S54) 


where {w™} is a set of choices of the eigenvalues = ±1 for all 2-WSO’s IT™, if exist, contained in the region A. 
Here 72 labels the non-contractible 2-cycles contained in the region A along which the 2-WSO W™ is constructed. In 
the presence of non-contractible 2 -cycles, the boundary conditions are also divided into different different topological 
sector whose label is given by the set choices {^^ 2 } for the eigenvalues of all the 2-WSO’s contained in the 
region A. Within each topological sector, the boundary conditions are related to each other, rendering the Schmidt 
coefficient only a function of and the {rc™}- In the absence of non-contractible 2-cycles, the corresponding 

labels drop out automatically. Also, in all the relevant cases, the numbers of permissible boundary conditions satisfy 
that if the region A is a union of two disjoint regions (say A^ and A^), the number of permissible boundary conditions 
on dA is a product of the number of permissible boundary conditions on dA^ and that on dA^. With these general 
understanding of the Schmidt coefficients, we can derive the all the relevant LRMI. The results are summarized as 
follows: 


I{D^) = 0, 

I{D^ X Ti) = 0, 

iiD^xT^y) = - ^ ^~ log ^~ I ^ log ^ 

aG{e,m} 

I{D^ X {T^y/D^)) = 0, 

1 - ip„ 1 I 1 ipg 1 + 


I{D^xT!y,) = - Y. Y 

aG{e,m} -yG{xy,yz,zx} 

X {TY/D^)) = - Y E 

a^{e,m} 'y^{xy,yz,zx} 


1 - 1 ^^. 1 - I 1 + TOy, 1 + (iTfu 


(S55) 


The generically non-vanishing values of I{D^ x —I{D^ x {T^y/D'^)) signifies the existence of the 2-membrane 
condensation. In contrast, the fact that 0 = I{D'^) = I{D^ x T^) = x Yy^) — I{D^ x {Tly^/DA indicates the 

absence of the 0-membrane, I-membrane and 3-membrane condensations. 
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